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We investigate the one-loop effect of new cliarged scalar bosons on tlie Higgs potential at finite 
temperatures in the supersymmetric standard model with four Higgs doublet chiral superfields as 
well as a pair of charged singlet chiral superfields. In this model, the mass of the lightest Higgs 
boson h is determined only by the D-term in the Higgs potential at the tree-level, while the triple 
Higgs boson coupling for hhh can receive a significant radiative correction due to nondecoupling one- 
loop contributions of the additional charged scalar bosons. We find that the same nondecoupling 
mechanism can also contribute to realize stronger first order electroweak phase transition than that 
in the minimal supersymmetric standard model, which is definitely required for a successful scenario 
of electroweak baryogenesis. Therefore, this model can be a new candidate for a model in which the 
baryon asymmetry of the Universe is explained at the electroweak scale. 



PACS numbers: 14.80. Da, 12.60.Fr, 12.60.Jv 

I. INTRODUCTION 

It has been confirmed that the numbers of matter and 
anti-matter in our Universe are asymmetric. The recent 
observation suggests that the ratio for baryon-to-photon 
is given by nb/n^ ~ (5.1 -6.5) x lO^io at the 95 % CL Q, 
where rif, is the difference in number density between 
baryons and anti-baryons and is the number density 
of photons. Understanding the mechanism of generating 
the baryon asymmetry of the Universe has been one of 
the most important issues in particle physics and cos- 
mology for a long time. In order to generate the baryon 
asymmetry from the baryon symmetric initial state, the 
following Sakharov's three conditions have to be satis- 
fied 0: 1) Baryon number nonconservation, 2) C and CP 
violation, 3) Departure from the thermal equilibrium. 

The electroweak gauge theory can naturally satisfy the 
above three conditions. The baryon number noncon- 
servation is realized by the sphaleron process at high 
temperature. The C violation can naturally occur in 
chiral gauge theories, and CP violating phases can re- 
main in the matter sector after rephasing. Finally, the 
circumstance of thermal non-equilibrium can appear if 
the electroweak phase transition (EWPT) is of strongly 
first order. This scenario is often called the electroweak 
baryogenesis 0, 0] ■ The most attractive feature of this 
scenario would be its testability. The scenario is neces- 
sarily related to the dynamics of electroweak symmetry 
breaking, so that it can be directly tested by measuring 
the Higgs boson properties at collider experiments. Cur- 
rently, Higgs boson searches are underway at the Teva- 
tron and the LHC. We expect that the physics of elec- 
troweak symmetry breaking will be clarified in near fu- 
ture. Then the models of electroweak baryogenesis can 
be experimentally tested. 



In the standard model (SM), it is found that the CP 
violation by the Cabibbo-Kobayashi-Maskawa matrix is 
quantitatively insufficient In addition, the require- 
ment of sufficiently strong first order EWPT predicts the 
Higgs boson mass to be too small to satisfy the constraint 
from the LEP direct search results @. Therefore, in or- 
der to realize the baryon asymmetry of the Universe we 
need to consider an extension of the SM. 

One of the viable models for successful electroweak 
baryogenesis would be the two Higgs doublet model 
(THDM) 0. The extension of the Higgs sector can in- 
troduce additional CP violating phases, and the quantum 
effect of extra scalar bosons in the Higgs potential makes 
it possible to realize sufficiently strong first order EWPT 
without contradictingthe LEP data for the mass of the 
Higgs boson. In Rcf. Q , the connection between the first 
order EWPT and the triple coupling for the lightest SM- 
like Higgs boson h (the hhh coupling) has been clarified. 
In the model with sufficiently strong first order EWPT, 
the hhh coupling constant significantly deviates from the 
SM prediction due to the same nondecoupling quantum 
effects of additional scalar bosons which make the first 
order EWPT strong. Such nondecoupling effects on the 
hhh coupling constant have been studied in Ref. . The 
scenario of electroweak baryogenesis by nondecoupling 
loop effects of extra bosons has also been applied in a 
TeV scale model where tiny neutrino masses, dark mat- 
ter and the baryon asymmetry of the Universe may be 
simultaneously explained [loj . 

It would also be attractive to consider the scenario of 
electroweak baryogenesis in the model based on super- 
symmetry (SUSY). SUSY is a good candidate of new 
physics, which eliminates the quadratic divergence in the 
one-loop calculation of the Higgs boson mass. The light- 
est SUSY partner particle in SUSY models with the R 
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parity can naturally be a candidate for the cold dark 
matter. 

In the minimal supersymmetric SM (MSSM), there 
are many studies to realize the electroweak baryogcnc- 
sis |ll| - [l3| . Currently, this scenario is highly constrained 
by the experimental data, especially the LEP Higgs mass 
bounds, leading to the tension between the lightest Higgs 
boson mass and the strength of the first order EWPT. 
Nevertheless, it is still viable for some specific mass spec- 
trum. According to Ref. [l^l, the strong first order 
EWPT is possible if < 127 GeV and m^^ < 120 GcV, 

where h is the lightest Higgs boson and ti is the lightest 
stop. To satisfy the LEP bound on rrih, the soft SUSY 
breaking mass for the left-handed stop should be greater 
than 6.5 TeV. The most striking feature of this scenario is 
that the electroweak vacuum is metastable and the global 
minimum is a charge-color-breaking vacuum, where the 
lifetime of the electroweak vacuum is found to be longer 
than the age of the Universe. 

The aforementioned tension in the MSSM baryogenc- 
sis can be relaxed by extending the Higgs sector. One of 
the possibilities is to add a gauge singlet field into the 
MSSM. So far, many studies on the electroweak baryo- 
genesis have been done in such singlet-extended MSSMs; 
i.e., the Next-to-MSSM [li], the nearly MSSM or the 
minimal non-MSSM the f/(l)'-extended MSSM [Hi, 
the secluded C/(l)'-cxtendcd MSSM [13, and so on. 
In this class of the models, the strong first order EWPT 
can be induced by the trilinear mixing term of the doublet 
and the singlet fields appearing in the tree-level Higgs po- 
tential. If this is the case, the mass constraints on the 
lightest Higgs boson and the light stop would be allevi- 
ated significantly. In the singlet-extended MSSM, how- 
ever, the vacuum structure is inevitably more compli- 
cated than the MSSM, giving rise to the unrealistic vacua 
in the large portion of the parameter space , espec ially 
electroweak baryogenesis-motivated scenario [ij, • 

In this Letter, we consider how the electroweak phase 
transition can be of sufficiently first order due to the 
nondecoupling effect of additional scalar bosons in an- 
other extended SUSY standard model, where a pair of 
extra doublet chiral superfields H3 {Y = —1/2) and H4 
{Y = +1/2) and a pair of charged singlet chiral super- 
fields Oi (y = -1-1) and 5I2 {Y = —1) are introduced 
in addition to the MSSM content [l^. Since the sin- 
glet fields are charged rather than neutral, they do not 
have the vacuum expectation values (VEVs) as long as 
the U{1)-EM is preserved. Therefore, at least at the tree 
level, there is no vacuum instability caused by the singlet 
field as opposed to the neutral singlet-extended MSSM 
as mentioned above. On top of this, this model may be 
motivated by SUSY extensions of the models with ad- 
ditional charged singlet fields in which neutrino masses 
are generated by radiative corrections [l^, HO]- In the 
present model, apart from the model in which an ad- 
ditional neutral singlet chiral supcrficld is added to the 
MSSM, there is no tree-level F-term contribution to the 
mass of the lightest Higgs boson, but there can be large 



H2 



Hz 
Hi 

^2 



Spin 



Spin 1/2 



11 = 



dR 



$1 = 
$2 = 



*1L = 



$4 



"I'sr, = 



-l-i 



SU(3) SU(2) U(l)., 



+ 1 



+ 1 



_ 1 
2 

+1 



2 

1 +1 
1 -1 



+ 



+ 



TABLE I: Property of chiral superfields (and their component 
fields) under the symmetries of the model. 



one-loop corrections to the triple Higgs boson coupling 
due to the additional bosonic loop contribution [l^ . We 
note that in the SUSY Higgs sector with four doublet 
chiral superfields (2]| . all the quartic interactions in the 
Higgs potential come from the D-term so that there are 
no large nondecoupling quantum effects on the hhh cou- 
pling [il. Similarly to the case of the non-SUSY THDM, 
these nondecoupling bosonic loop contributions can also 
make first order phase transition stronger. We here show 
that the EWPT can be of sufficiently strong first order 
in this model. 



II. MODEL 

We consider the model with the chiral superfields 
shown in Table 1 [l^. The symmetries of the model 
are standard gauge symmetries. In addition, we impose 
a discrete Z2 symmetry for simplicity. Although the Z2 
symmetry is not essential for our discussion, the symme- 
try works for avoidin g th e fiavor changing neutral current 
at the tree level [2l| - |24| . Furthermore, we assume that 
there is the R parity in our model. 

The superpotential is given by 

W = {yuTUtH^ ■ Q, + {yarD'TH, ■ Q, + [y^f'ElH^ ■ L, 

- ^iH, ■ H2 - ^i' H3 ■ Hi - fianiQ2- (1) 
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The soft-SUSY-breaking terms are given by 
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Csott = --{MiBB + A'hWW + MiGG) 

+ (Ae)'^"g|,,<I>i • 

• $3 + (A2)a;2"$2 • $4 + h-c.} . (2) 
From W and -Cg^j^, the Lagrangian is constructed as 

^ '^kinetic ~t~ -^gaugc— matter 



l_^W_ 

2 dipidtpj 

2 1 



'(pLi ■ i'Lj + /i-C. 



dW 



-'Soft 



(3) 



where tpi and tpLi are respectively scalar and fermion 
components of chiral superfields, and T^^ and ga rep- 
resent generator matrices for the gauge symmetries and 
corresponding gauge coupling constants. 

The scalar component doublet fields are parameter- 
ized as 
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$2,4 



4>tA 

^(V32,4 + /l2,4 + 



ia2,4) 
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where (pi are classical expectation values, hi are CP-even, 
tti are CP-odd and (j>^ are charged scalar states, where ipi 
are classical expectation values, hi are CP-even, are 
CP-odd and are charged scalar states. We use the 
effective potential method to explore the Higgs sector. 
At the tree level, the effective potential for the Higgs 
fields is given by 

K)(<^1,<^2,<^3,<^4) 
= o'^a'^a + T^iB^Lpiipo^ + B'fl'Lp3Lpi + h.C.) 



32 



/ 2 2,2 2\2 



(5) 



Using the effective potential, the vacuum is determined 
by the stationary condition as 



dipi 



= 0. 



(6) 



We assume that the Z2 odd Higgs bosons do not 
have the VEVs (^3 = V4 = 0) for simplicity, and 
we set \/ vl + V2 = f (~ 246 GeV) and introduce 
tan /? = V2/V1. At the tree level, W3 = 1)4 = is 
guaranteed by requiring the nonnegativc eigenvalues of 
{d'^Vo/dtpidLpj)^^ .^=o = 3,4), i.e.. 



s'V^ > 0, 



> 0. 



(7) 



In the following, we exclusively focus on the ((^1,(^2) 
space. 

For the Z2 even scalar states, after the symmetry 
breaking we have five physical states as in the MSSM; 
i.e., two CP-even h and H, a CP-odd A and a pair of 
charged scalar bosons. The tree level mass formulae 
for these scalar states coincide with those in the MSSM. 

The mass cigenstates for the Z2 odd charged scalar 
states (j)'-^ and (f)'^ are obtained by diagonalizing the com- 
ponent fields of doublet scalar fields $3 and $4, and fjf 
and from the charged singlet scalar fields uji and UJ2- 
Their field dependent masses are given by 



<Pl,2 


1 

~ 2 


ml 




1 

~ 2 





ml + ml + - (lAilV? + IA2IV2) T 



^- + 2(I-^i|V? + |A2|V2)t/D^ 

(8) 



where 
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(10) 



Mass parameters in the above formulas are 4 = 
^^3.4 + ImT and rril^_ = _ + \^^\2^ ^he physi- 
cal masses can be obtained by replacing tpi by Vi, the 
values at the vacuum. On the other hand, masses of 
the additional Z2 odd neutral scalar bosons do not re- 
ceive the contributions from the |Ai|^-terms, and only re- 
ceive those from the D-tcrm and the soft-SUSY-breaking 
terms. Their loop effect on the effective potential are 
small, so that we neglect them in our later calculations. 

The field dependent masses of the Z2 odd charginos 
are given by 



m'i,± 
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(11) 



4 



where 
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{\Xl\'\Vl\'-\X2f\^2\'] 



(12) 



The physical masses can be obtained by replacing ipi by 
Vi. Notice that the masses do not vanish even when the 
invariant mass parameters /x' and /io are taken to be zero. 

Since it is known that radiative corrections on the 
Higgs sector arc very important to study the EWPT, 
we here focus on the one-loop contribution. The vacuum 
at the one-loop level is also determined from Eq. © with 
the one-loop corrected effective potential. The one-loop 
correction to the effective potential at zero temperature 
is given by 
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FIG. 1: The masses of the Z2-0AA scalars and Z2-oAd 
charginos as a function of A2. The masses of heavy Z2-even 
neutral scalars H and A are also plotted. Model parameters 
are taken as them shown in Eq. p7p . The Xt dependence of 
these mass spectrum is negligible. 



where Vi is regularized in the DR-scheme, a is the de- 
grees of freedom of the species i, M is a renormalization 
scale which will be set on m^°'''. 

For the zero temperature T = 0, the one-loop corrected 
mass matrix for the CP even neutral bosons can be calcu- 
lated from the effective potential. We here consider the 
simple case such that B' = Bq = = /.tfj = in order 
to switch off the mixing effects. By using the effective 
potential method, the MSSM-like CP even Higgs boson 
mass matrix (M^)ij with the leading Xf 2 contributions 
is given as 
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FIG. 2: The Z2-even lightest Higgs boson mass as a function 
of A2. From the top to the bottom, Xt/Mq = 2.0, 1.2 and 0.6. 



where 1 , 2 are labeled as 

n2 



Xi 



X2 



nir 



TTlr 



C2P 
C2l3 



X2„2„2 



-^2^13 2 
-V 



AlC^^2 



-2 , / 2 2 ^ , '^Z*/? 2 

=m_ + [mz - m^)c2fj H —v 



-2 I 2 2 \ I 

=771+ - (m^ - m,y)c2/3 + 



Xlc' 



%2 



(15) 



The rcnormalized mass of the lightest Higgs boson h is 



then calculated for to^ 3> mz as 
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at the leading X\ 2 contributions, where the one-loop con- 
tribution in the MSSM is mainly from the top and stop 
loop diagram [2^. 

Now we quantify the magnitude of the radiative cor- 
rections of the Z2-odd particles on m/j. The input pa- 
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rameters are fixed as follows. 

Tree : tan/3 = 3, m^i = 500 GcV; 

1-loop (MSSM) : Mq = = Mf = 1000 GeV, 
/i = Ah = 2Mi = 200 GeV, 
At = At = Xt + /i/ tan/3; 

1-loop ($'±2. -^1 2, /i' = = Bn = S' = 0, 
rn:^ = rn^ = (500 GeV)^ 
m^_^ml^ (50 GeV)^ (17) 

We note that m^/± < m^;± and m^^± < m^^± in this case. 
The mass spectrum of extra Z2 odd charged scalars and 
charginos is displayed in Fig. [1] On this parameter set, 
'm'i_i+, TO?,+ and m^,+ get a significant contribution from 

A2 . Then their masses become larger for the greater value 
of A2 . Since the mass parameters m\ and ifi^_ are taken 
to be small, large mass values of m^'± and TOq± yield the 
large nondecoupling effects which can make the EWPT 
strongly first order. 

Fig.[2]shows the predicted value of mh as a function of 
A2 varying Xt/Mq = 2.0, 1.2 and 0.6 from the top to the 
bottom. We can see that rrih monotonically decreases as 
A2 increases, which is in contrast with the top/stop loop 
effects. 

The coupling constants Ai and A2 are free parame- 
ters of the model. Its magnitude, however, is bounded 
from above by the condition that there is no Landau 
pole below the given cutoff scale A. As we are interested 
in the model where the first order EWPT is sufficiently 
strong, we allow rather larger values for these coupling 
constants, and do not require that the model holds un- 
til the grand unification scale. A simple renormaliza- 
tion group equation analysis tells us that for assuming 
A = 2 TeV, 10 TeV or 10^ TeV, the coupling constant 
can be taken to be at most A2 ~ 2.5, 2.0 or 1.5, re- 
spectively. Above the cutoff scale A, the model may be 
replaced by a strongly coupled supersymmetric theory 
with UV completion as described by the scenario such as 
in the fat Higgs model [26j . 

III. ELECTROWEAK PHASE TRANSITION 

The nonzero temperature effective potential is 

V,{v,,^2-T)=Y^c~Ib,f (^) , (18) 

where B{F) refer to boson (fermion) and Ib,f take the 
form 

lB,Fia^)--J dx ln(lTe-^"'+"J. (19) 

Since the minimum search using Ib,f is rather time- 
consuming, we will alternatively use the fitting functions 



of them that are employed in Ref. [T^. More explicitly, 

N 

/B,^(a2) = e--^c^./a", (20) 

n=0 

are used, where c^^ are determined by the least square 
method. For N = 40, \lB,F{a'^) - lB,F{a'^)\ < lO'^ for 
any a, which is sufficient in our investigation. Since the 
nonzero modes of the thermal effective potential give 
corrections to the 2-point self energy at high tempera- 
tures, we will resum them to make the perturbative anal- 
ysis more reliable p7| . 

For an electroweak baryogenesis scenario to be suc- 
cessful, the sphaleron rate in the broken phase should be 
smaller than the Hubble constant. Conventionally, this 
condition is translated into 

where Tq is the critical temperature, vc is the Higgs 
VEV at Tc, and C is a 0(1) parameter. To obtain C 
within a better accuracy, the sphaleron energy and zero- 
mode factors of the fluctuations around the sphaleron 
must be evaluated. In the SM, the sphaleron energy is 
simply a function of the Higgs boson mass. As the Higgs 
boson becomes heavier, the sphaleron energy gets larger 
as well [2^, leading to the smaller ( [2^. In this model, 
on the other hand, ^ depends on more parameters. For 
simplicity, we here take C = 1; which is often adopted as 
a rough criterion in the literature. 

In our analysis, Tc is defined as the temperature at 
which the effective potential has the two degenerate min- 
ima. We search for Tc by minimizing 

(22) 

where the field-dependent masses are modified by adding 
thermal corrections, and Eq. ([20]) is used in Vi{ipi,ip2',T). 

In Fig. [3] Vc and Tc are plotted as a function of A2 
in the fight h scenario {Xt/Mq = 0.6): see Fig. H The 
sphaleron decoupling condition (|2ip can be fulfilled for 
A2 > 1.6 due to the nondecoupling effects coming from 
(j)'^ and VL^. 

We also evaluate vc and Tc in the heavy h scenario 
{Xt/Mq = 2.0) as shown in Fig. SI The sphaleron de- 
coupling condition can be satisfied for A2 > 1.8. Though 
the parameter region is a bit narrower than the light 
Higgs scenario, the lightest Higgs boson mass as large as 
134 GeV is still consistent with the decoupling condition. 



IV. DISCUSSIONS 

The several comments on the current analysis are in 
order. 
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FIG. 3: vc and Tc vs. X2 with Xt/Mg = 0.6. The other 
input parameters are the same as in the Fig. [2] The sphaleron 
decouphng condition (|2ip can be satisfied for A2 > 1.6. 




FIG. 4: vc and Tc vs. A2 with Xt/Mg = 2.0. The other 
input parameters are the same as in the Fig. [21 The sphaleron 
decoupling condition (|2ip can be satisfied for A2 > 1.8. 



• According to a study of the sphaleron decoup ling 
condition in the MSSM, it is found that ( ~ 1.4 [l3| 
which is 40% stronger than one we impose in our 
analysis. The similar value may be obtained in this 
model as well. It should be emphasized, however, 
that even if wc take ^ = 1.4 for the sphaleron de- 
coupling condition, a feasible region exists for the 
relatively large A2, for example, A2 > 2.2 even in 
the heavy h case. The cutoff scale A is rather low in 
this case but still it is around the multi-TeV scale. 

• In this model, the light stop scenario is one of the 
options for the successful electroweak baryogenesis. 
Same as the scenario in the MSSM, the strength 
of the first order EWPT can get enhanced if the 
(almost) right-handed stop is lighter than the top 
quark, enlarging the possible region. 



• Similar to the usual MSSM baryogenesis scenario, 
the charginos and/or the neutralinos can play an 
essential role in generating the CP violating sources 
as needed for the bias of the chiral charge densities 
around the Higgs bubble walls. In addition to this, 
the Z2 odd charginos Xi^2 ^^^^ do the job for 
the successful baryogenesis. 

• Since Ai^2 and Ai,2Mn small in our parameter 
choices, the charge breaking does not occur at the 
tree level. In addition, the Z2 symmetry is not 
broken spontaneously at the tree level, because m|, 
ml and B' fi' we take here satisfy Eq. ([7|). The 
potential analysis beyond the tree level is out of 
scope in this Letter. It will be our future task. 

• If the Z2 symmetry is exact and unbroken after 
the electroweak symmetry breaking, the lightest 
Z2 odd particle in our model can be a candidate 
of cold dark matter if it is electrically neutral, in 
addition to the lightest supersymmetric particle. If 
one of the extra neutral scalar bosons is the lightest 
Z2 field, its phenomenological property and exper- 
imental constraints would be similar to those for 
the supersymmetric extension of the inert doublet 
model [13] ■ A neutralino from the extra doublets 
may also be a candidate for dark matter. 

• Finally, wc comment on the phenomenological pre- 
dictions of this model. First of all, the nonde- 
coupling effect of the extra Z2 odd charged scalar 
bosons on the finite temperature effective potential 
is an essentially important feature of our scenario 
in order to realize strong first order phase tran- 
sition. The same physics affects the triple Higgs 
boson coupling with a large deviation from the SM 
(MSSM) prediction as discussed in Ref. Q, Such 
deviation in the triple Higgs boson coupling can be 
15-70 % H, [1^, and we expect that they can be 
measured at the future linear collider such as the 
ILC or the CLIC. Second, in our model, in order to 
realize the nondecoupling effect large, the invariant 
parameters fi' and /lo are taken to be small. Conse- 
quently, the masses of extra charginos are relatively 
as hght as 100-300 GcV. 



V. CONCLUSIONS 

We have discussed the one-loop effect of new charged 
scalar bosons on the Higgs potential at finite tempera- 
tures in the supersymmetric standard model with four 
Higgs doublet chiral superfields as well as a pair of 
charged singlet chiral superfields. We have found that 
the nondecoupling loop effects of additional charged 
scalar bosons can make first order EWPT strong enough 
to realize successful electroweak baryogenesis. We, 
therefore, conclude that this model can be a new good 
candidate for a successful model where the baryon 
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asymmetry of the Universe is explained at the elec- 
troweak scale. The detailed analysis for the collider 
phenomenology will be shown elsewhere. 
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